Comment on "Nonuniversality of Transverse Momentum Dependent Parton 

Distributions at Small x" 
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We point out that the analysis in arXiv: 1003.0482 actually verifies the universality of transverse 
momentum dependent quark distributions at small x, which supports the observation in our earlier 
work arXiv:0904.4150 Once the gluon exchanges responsible for the additional infrared divergences 
are factorized and summed into an exponential factor, it can be expressed as a matrix element 
of Wilson line operators, and treated as a nonperturbative input to the kr factorization of pA 
collisions. The same exponential factor has been extracted from pp collisions in arXiv:0904.4150 , 
rendering possible experimental constraints on its behavior from pp processes, and predictions for 
pA processes. 



In a recent letter Xiao and Yuan analyzed the ini- 
tial and final state interaction effects in dijet-correlations 
in pA collisions at small x, and demonstrated the nonuni- 
versality of the transverse momentum dependent (TMD) 
quark distributions in a scalar-QED model. The nonuni- 
versality is caused by the infrared divergence attributed 
to vanishing invariant mass squared — k\ of a Glauber 
gluon, kj_ being the gluon transverse momentum. As a 
comparison, the all-order amplitudes for the TMD quark 
distributions (see Eq. (14) and Eq. (39) in @) 
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were extracted from the deep inelastic lepton-nucleus 
scattering (DIS) and the dijet-correlation process p+ A — > 
Jetl + Jet2 + X, respectively. The factor G(R X + rj_) — 
G(Rl) with G(Rl) = Ko(XR±)/(2tt), K being the 
Bessel function, and A an infrared regulator, approaches 
-W(r x ,R±) = (l/27r)ln(i?_ L /|R L +r±\) in the A -> 
limit. The additional exponential factor e~' l992G ( R±S) 
arises from the summation of the Glauber divergences to 
all orders. The TMD quark distribution for the nucleus 
A from the dijet correlation was then given by 
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with = R± + r± — R' ± . It was claimed that the 
Glauber factor breaks the universality of the TMD quark 
distribution (l|. 

We first point out that the same Glauber factor e~ lS ^ 
has been derived in our earlier work on the di-hadron 
production in pp collisions at small x in the same model 
[3j, with the exponent S from the one-loop correction, 
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Equation (O leads to G(Rj_) in [l[ with the impact pa- 
rameter b. The difference is that we did not differentiate 
the coupling constants gi and g associated with the par- 
ton lines from the projectile and from the target, respec- 
tively. The Feynman rule in Eq. ([3]) was constructed by 
means of the eikonal approximation for the Glauber glu- 
ons at small x Hence, it is straightforward to obtain 
the definition of the Glauber factor 
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where the Wilson line operator is written as 
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The expression of W + is similar, but with the light-cone 
vector = (0,1,0) being replaced by = (1,0,0), 
and 32 by g. The net effect of W_(Rj_; — oo) t T4 / -(Rj_; oo) 
and W+(0; oo)W+(0; — oo)^ demands the vanishing of the 
components k + and k~ of the gluon momentum, respec- 
tively. It is easy to show, by expanding the Wilson line 
operators order by order, that Eq. ^ reproduces Eq. ([3]). 

Obviously, Eqs. (fTJ) and ([2]) indicate that the Glauber 
divergences are factorizable in the impact parameter 
space. Once the factorization is achieved, the Glauber 
factor can be treated as an additional nonperturbative 
input to pA collisions, just like the soft factor introduced 
in the kr factorization 0] ■ The same Glauber factor was 
derived in Q and [l], , rendering possible experimental 
constraints on its behavior from pp processes, and pre- 
dictions for pA processes. Based on the observations that 
the Glauber factor can be expressed as a matrix element 
of the Wilson line operators, treated as a convolution 
piece in the kx factorization, and constrained from some 
processes, we postulate that the analysis in [l[ actually 
verifies the universality of TMD quark distributions at 
small x as claimed in [3J. 
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